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Abstract 

Chromatin remodeling machineries are abundant and 
diverse in eukaryotic cells. They have been involved 
in a variety of situations such as histone exchange and 
DNA repair, but their importance in gene expression 
remains unclear. Although the influence of nucleosome 
position on the regulation of gene expression is gener¬ 
ally envisioned under the quasi-equilibrium perspective, 
it is proposed that given the ATP-dependence of chro¬ 
matin remodeling enzymes, certain mechanisms neces¬ 
sitate non-equilibrium treatments. Examination of the 
celebrated chromatin remodeling system of the mouse 
mammary tumor virus, in which the binding of tran¬ 
scription factors opens the way to other ones, reveals 
that breaking equilibrium offers a subtle mode of tran¬ 
scription factor cooperativity, avoids molecular trapping 
phenomena and allows to reconcile previously conflicting 
experimental data. This mechanism provides a control 
lever of promoter responsiveness to transcription factor 
combinations, challenging the classical view of the uni¬ 
lateral influence of pioneer on secondary transcription 
factors. 

Keywords Chromatin remodeling; cooperativity; 
transcription factor; MMTV; glucocorticoid receptor. 

1 Introduction 

The importance of ATP-dependent machineries remod¬ 
eling chromatin by actively moving nucleosomes rela¬ 
tively to DNA, remains puzzling. Beside their possi¬ 
ble structural roles in chromatin organization, nucleo¬ 
some repositioning, histone exchange and DNA repair, 
a role in transcriptional cooperativity is proposed here. 
The specificity and intensity of gene expression is gov¬ 
erned by interactions between regulatory DNA sequences 
(cis-regulators) and various trans-acting factors (tran¬ 
scription factor proteins (TFs) and non-coding RNAs). 
The occupation of a gene promoter by these trans¬ 
regulators involves both micro-reversible and micro- 
irreversible steps. Micro-reversible binding processes 
can lead to sigmoidal concentration-dependent response 


through classical multimeric cooperativity (Bolouri and 
Davidson, 2002; Michel, 2010). The role of nucleosomes 
has also been examined from the micro-reversible per¬ 
spective (Dodd et ah, 2007; Segal and Widom, 2009; 
Mirny, 2010). The rapid equilibration of these thermally- 
driven phenomena, relatively to the slow changes of cel¬ 
lular components, simplifies the definition of the input 
functions used in gene network modeling (Bintu et ah, 
2005; Michel, 2010). But promoter occupancy also in¬ 
volves some micro-irreversible transitions such as chro¬ 
matin remodeling and active dissociation processes. Pre¬ 
cisely, it is proposed in the present study that insert¬ 
ing micro-irreversible steps in the process of promoter 
saturation, offers additional possibilities of potent coop¬ 
erativity. A single example has been selected because 
it remarkably illustrates how micro-irreversible transi¬ 
tions can generate a refined discernment in gene ex¬ 
pression. In this example, the micro-irreversible step 
corresponds to the energy-dependent phenomenon of 
chromatin-remodeling, in which the position of DNA 
around nucleosomes is modified. By this way, upon bind¬ 
ing to DNA, a first TF directs the accessibility to DNA 
of other ones. This mechanical activity allows: (i) coop¬ 
erativity between non-physically interacting TFs and (ii) 
constitutively expressed TFs to participate to conditional 
induction. One of the most celebrated chromatin remod¬ 
eling system is provided by the thoroughly documented 
Mouse Mammary Tumor Virus (MMTV) promoter. 

2 Breaking hierarchical polymer¬ 
ization is necessary to maintain 
molecular dynamics 

The assembly of macromolecular complexes generally 
proceeds in a hierarchical manner in the cell. For ex¬ 
ample, a component C which cannot bind to the isolated 
components A and B, can bind to a pre-associated com¬ 
plex AB. Hierarchical binding chains such as A -|- H ^ 
AB, +C ^ ABC, -\-D ^ ABCD etc, are often involved 
in the building of multi-molecular complexes, but are less 
compatible with the dynamic and reactive behaviours of 
solubles components. Indeed, in equilibrium conditions, 
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the chain written above leads to the trapping of the early 
components in the complexes as long as the late compo¬ 
nents are present. This phenomenon can exist for TF 
binding to gene promoters. It can hold for example, for 
the successive binding steps observed in equilibrium con¬ 
ditions between the TodT TFs and the series of TodT- 
binding sites juxtaposed in the Tod gene promoter, that 
has been proposed to be mediated by DNA conformation 
changes (Lacal et ah, 2008). Beside this puzzling situa¬ 
tion of equilibrium allostery, the hierarchical binding of 
transcription factors in equilibrium conditions is also pos¬ 
sible in the case of the large eukaryotic preinitiation com¬ 
plexes made of the so-called general transcription factors 
GTFs (Michel, 2010). But hierarchical relationships have 
also been reported for non-interacting isolated TFs in ab¬ 
sence of any trapping phenomenon. To allow reconciling 
hierarchical binding and absence of trapping, one should 
postulate the possibility to break equilibrium. This sit¬ 
uation is well illustrated by the case of the occupation 
of the MMTV promoter involving micro-irreversible pro¬ 
cesses, thoroughly documented but not yet clearly un¬ 
derstood. In section 4, this system will be analysed un¬ 
der the classical equilibrium assumption. Then, in sec¬ 
tion 5, a non-equilibrium scheme will be proposed based 
on hypotheses built from MMTV experimental data, in 
which the equilibrium-breaking machines are the ATP- 
dependent SWI/SNF chromatin remodeling enzymes. 

3 Data obtained for the MMTV 
promoter occupancy are ir¬ 
reconcilable from the time- 
reversible perspective 

A central piece of data about MMTV expression is the 
role of nucleosomes in the mutual influence between 
the glucocorticoid receptor (GR) and a group of TFs 
(NFl/Oct). Although the activation of MMTV by GR 
and NFI/Oct-1 seemed clear in the initial reports, dis¬ 
crepancies appeared in following studies. The basis 
of glucocorticoid hormone-induced MMTV regulation is 
that GR has an initiating role, triggered upon hormone 
binding (stress hormone corticol or corticosterone) and 
subsequently amplified by NFI and Oct-1. This sequen¬ 
tial action is dependent on the position of nucleosomes on 
DNA, since it is not observed with naked DNA (Richard- 
Foy and Hager, 1987; Archer et ah, 1992; Ghavez and 
Beato, 1997). The repositioning of nucleosomes triggered 
by GR, probed by nuclease or chemical mapping, leads 
to the exposure of the NFI and Oct-1 binding sites and 
is mediated by SWI-SNF ATPases (Fryer and Archer, 
1998). The different roles of GR and NFI in initiating 
and amplifying transcription respectively, are explained 
by their differential mode of interaction with chromatin: 
GR can bind to DNA wrapped around nucleosomes, con¬ 


trary to NFI which requires a fully accessible double he¬ 
lix (Eisfeld et ah, 1997). This DNA-binding hierarchy, 
first of hormone-bound GR and then of NFl/Oct, of¬ 
fers a powerful opportunity of cooperativity. Indeed, GR 
is inducible but not very potent contrary to the couple 
NFl/Oct. As NFl/Oct can access DNA only upon bind¬ 
ing of GR, the promoter activity can become strongly 
sigmoidal, particularly if NFl/Oct is transcriptionally 
more potent than GR. Sigmoidal responses are gener¬ 
ally due to decreased responsiveness to low signals. This 
is the case for the MMTV promoter in which TFs are 
prevented to bind at low concentration. But this ele¬ 
gant mechanism has then been clouded in the following 
reports, which introduced new actors and revised the hi¬ 
erarchy of binding of GR and NFI. In sharp contrast 
with the earlier articles, NFI and Oct-1 binding sites 
have been shown to preset chromatin prior to GR binding 
(Belikov et ah, 2004). The picture is thus more complex 
than supposed previously and the mutual influence be¬ 
tween GR and NFI for binding to the MMTV promoter 
is now described as dualistic (Hebbar and Archer, 2007), 
blurring the logic of this system. The same apparent 
paradox has been pointed for the relationships between 
purported ” pioneer” and secondary transcription factors 
(Gaizzi et ah, 2014). In fact, the strict dependence on 
the previous fixation of a factor to allow the fixation of 
another factor can be alleviated if introducing an ad¬ 
ditional step of chromatin remodeling (Fig.l). It will 
be shown that the cooperative relationships between GR 
and NFI, which are unclear when examined only from a 
time-reversible perspective, can be usefully reconsidered 
from a non-equilibrium perspective (section 5), but the 
outcomes of equilibrium modeling is first examined below 
for comparison. 

4 Equilibrium modeling of hierar¬ 
chical MMTV promoter occu¬ 
pancy 

In the simplest hierarchical modeling scheme assum¬ 
ing micro-reversibility (Fig.la), MMTV transcription is 
stimulated by two groups of transcription factors GR 
(named A) and NFl/Oct (named B) (Fig.l). A and B 
bind to the MMTV promoter (P) through their DNA- 
binding domain (DBD) in a hierarchical manner, but 
once bound to DNA, they are supposed to stimulate 
transcription in an independent and additive manner, 
through their activation domain (TAD). In these con¬ 
ditions, the fractional activity (F) ranging from 0 to 1, 
of the MMTV promoter, is described in Eq.(l). 

p(A) kA+p{B) ks 

In this equation, and ks are the maximal frequen¬ 
cies at which A and B, when bound to DNA, recruit tran- 
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scription machineries, thereby initiating multiple rounds 
of transcription. These frequencies should be weighted 
by the probabilities of presence of A and B on the pro¬ 
moter, written p{A) and p{B) respectively (with small 
letters p to not be confused with the promoter P). 
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The probabilities p{A) and p{B), equivalent to frac¬ 
tional occupation times, can be formulated through an 
Adair approach, as the ratio of occupied over total bind¬ 
ing sites, which can be expressed as concentrations in er- 
godic conditions, by enumerating the possible promoter 
states. 


^ [PA] + [PAB] 

’ [Po] + [PA] + [PAB] 


(2a) 


and, given that B is supposed to access DNA only when 
A is present. 


p{B) 


[PAB] 

[Po] + [PA] + [PAB] 


(2b) 


Using Po as a reference, Eq.(2) can be converted into 
p{A) = Ka[A]{1 + Kb[B])/D (3a) 


p{B)=Ka[A]Kb[B]/D (3b) 

with 

D = 1 + Ka[A]+Ka [A] Kb [B] (3c) 


In this scheme, GR is prevented to dissociate from a 
saturated promoter. This trapping effect which can ap¬ 
pear puzzling, is inherent to the equilibrium modeling of 
sequential cooperativity, but such a trapping of GR on 
the MMTV promoter is not consistent with the observa¬ 
tion that GR can escape DNA whatever the chromatin 
configuration (Fletcher et ah, 2000; Hager et ah, 2000). 
To avoid this problem, one can imagine an alternative 
scenario (Fig.lb), in which GR dissociation does not re¬ 
quire the absence of NFl. Two different equilibrium con¬ 
stants are defined for GR {Ka and K'^^ = l/AT^^), to take 
into account the different chromatin states. The occupa¬ 
tion probabilities of the A and B sites are respectively: 

p{A) = Ka[A]{1 + Kb[B])/D (4a) 

p{B) = Ka[A]Kb[B]{1 + K',J[A])/D (4b) 

with 

D = l + Ka[A]{1 + Kb[B]{1 + K'^A/[m (4c) 


Figure 1. Different models to explain the hierarchical occu¬ 
pation of the MMTV promoter (P) by GR (named A) and 
NFl (named B). The schemes (a) and (b), comply with the 
principle of microscopic reversibility but not the schemes (c) 
and (d). In (c), ko is the rate of chromatin opening driven 
by SWI/SNF ATPases and kc is the rate of chromatin clos¬ 
ing, driven by stabilization of DNA bending. B cannot bind 
to P because of inappropriate nucleosome positioning, while 
A can bind to both P and P' with different constants. In 
the scheme (d), chromatin closing can occur only when the 
promoter is free of any TF. 


But the phenomenon of trapping now concerns NFl, 
possibly for long periods in case of removal of the glu¬ 
cocorticoid hormone. Though puzzling, this possibility 
would be consistent with the observation that NFl is 
present on DNA prior to hormone addition (Hebbar and 
Archer, 2003; Belikov et ah, 2004). But it remains to 
explain why the MMTV promoter would be inactive in 
spite of the continuous presence of NFl. A possible ex¬ 
planation could be that NFl doesn’t work as long as it 
is prevented to recycle on DNA, according to the model 
of one-shot TFs like ATF6 (Michel, 2010). This system 
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would still conform microscopic reversibility, but the no¬ 
tion of equilibrium would become shaky since it is no 
longer dynamic after disappearance of A. To propose 
a more plausible formulation of the MMTV promoter 
occupancy, relieved from any trapping effect, a micro- 
irreversible steps should be introduced in the system. 


5 Non-equilibrium modeling of 
the MMTV promoter occu¬ 
pancy 


To not recourse to trapping phenomena which are not 
experimentally verihed for MMTV, one should postulate 
another mode of cooperativity, liberated from the micro¬ 
reversibility constraints. Energy inputs obviously exist 
in the system and are provided by ATPases (SWI/SNF), 
recruited by DNA-bound GR (Fryer and Archer, 1998). 
Among the different micro-irreversible mechanisms that 
can be imagined, the model shown in Fig.lc is an at¬ 
tempt to reconcile the more experimental data as pos¬ 
sible, in a novel scheme as simple as possible. Chro¬ 
matin remodeling can be triggered and reversed dynami¬ 
cally, according to the well established reversibility (in its 
traditional acceptation) of hormone-induced nucleosome 
positioning (Belikov et ah, 2001) and to the dynamic in¬ 
teraction of remodeling complexes with the MMTV pro¬ 
moter (Johnson et ah, 2008). The spontaneous nucleo¬ 
some repositioning from P' to P is dictated by the intrin¬ 
sic bendability of DNA sequence and can be considered 
as nearly micro-irreversible. Since transient behaviours 
following GR activation can be neglected for the result¬ 
ing gene expression, a steady state treatment is sufficient 
for modeling this system. A reasonable additional hy¬ 
pothesis is that the time scales are different between the 
rapid equilibration of the TFs with the promoter, and 
the slower dynamics of micro-irreversible chromatin re¬ 
modeling, for opening site for B (ko) and for closing it 
(fee). The time scale separation hypothesis is not always 
applicable, but is justihed in the present case, given the 
rapid equilibration of the TFs with the MMTV DNA, 
suggested by the short turnovers of GR (12 milliseconds) 
evidenced by fluorescence recovery after photobleaching 
(FRAP)(Sprague et ah, 2004). This condition allows to 
use the approach of (Gha, 1968), mixing in the same 
treatment rate constants (time-dependent) and equilib¬ 
rium constants (time-independent). In this method, sev¬ 
eral groups of rapidly equilibrated species are defined 
using equilibrium constants. They correspond in the 
present case to the two chromatin states of the MMTV 
promoter, which will be written EP and EP' (Fig. lc,d). 
Two DNA-binding constants Ka and K'^ are postulated 
for A depending on the chromatin state, but K'^ is not 
affected by the presence or not of B, given that these 


TFs do not directly interact with each other. 

[EP] = [Po] + [PA] (5a) 

[EP'] = [P^j + [P'A] + [P'P] + [P'AP] (5b) 


with 


and 


[^o] 1 

[SP] Dp 

[PA] ^ Ka[A] 
[EP] Dp 


Dp = 1 + Ka[A] 


For the P' states: 


[P^] _ 1 

[EP'] Dp, 

[P'A] _ K'JA] 

[EP'] Dp, 

[P'P] _ Kb[B] 

[EP'] Dp, 

[P'AB] _ K'JA]Kb[B] 
[EP'] Dp, 

with 

Dp, = il + K'A[A])il + KB[B]) 


(6a) 

(6b) 

(6c) 

(7a) 

(7b) 

(7c) 

(7d) 

(7e) 


When gathering the promoter states, the probabilities 
p{A) and p(A) that P is occupied by A and P respec¬ 
tively, can be defined as follows 

p(A) = p(P n A) -I- p(P'n A) (8a) 

equivalent to 

p(A) = p(A|P) p(P) + p(A|P') p{P') (8b) 

and 

p{B) = p{P' n P) (given that p{P n P) = 0) (8c) 


where 


p{B) = p{B P') p{P') 


p{A\P) 


[PA] 

[EP] 


Ka[A] 

1 + Ka[A] 


p(A|P') 


[P'A] -f [P'AP] 

ppi 


K'aIA] 

1 + K'^[A] 


p{B]P') 


and 


[P'P] -f [P'AP] 

ppi 


Kb[B] 

1 + Kb[B] 


p{P) 


[P] 

[P] + [P'] 


and p{P') 


l-p(P) 


(8d) 

(8e) 

(8f) 

(8g) 

(8h) 
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[P] and [P'] are the amounts of time in which the pro¬ 
moter is in the P and P' states, which can be deduced 
from the steady state balance. If supposing, to agree 
with experimental observations, that the restoration of 
the basal chromatin can occur only when the promoter 
is not occupied by B (Fig.lc), then, the steady state can 
be written: 

ko[P] piA\P) = k,[P'] (1 - piB\P')) (9) 

which yields, using the values defined previously, 

^ k,il + KA[A]) 

[P'] koKA[A]{l + Kb[B]) 

leading to 

piA) = Ka[A] (^k, + koK'A[A] 
and 

p{B) = koKA[A]KB[B]/D 

where 

D = kcil + Ka[A]) + koKA[A]{l + Kb[B]) 


( 10 ) 

/D 

( ll a) 

( ll b) 

( ll c) 


These results are then incorporated in Eq.(l). The 
capacity of this system to generate sigmoidal curves is 
due to the products of the concentrations of A with it¬ 
self (in Eq.(lla)) and with B (in Eqs(lla) and (11b)) 
if assuming a double time scale separation: (i) between 
DNA/TF interactions and chromatin remodeling kinet¬ 
ics, as previously postulated and (ii) between chromatin 
remodeling and gene product concentration changes. 


6 TF concentration-dependence 
of the promoter activity 

The important parameters to evaluate are the sigmoidic- 
ity and sensitivity of the promoter activity to TF con¬ 
centration changes. Sigmoidicity is classically obtained 
when the TFs should multimerize for binding to DNA. 
This is the case for GR which is active as a dimer. But 
to focus on the specific source of cooperativity provided 
by the mechanism examined here, dimerisation will be 
ignored and the TFs will be considered as preformed 
dimers. For easier analyses, the fractional promoter ac¬ 
tivity equations will be adimensioned by setting some 
constants. The ratio of the transcriptional strength of 
DNA-bound B and A is 7 = kB/kA- In the mechanisms 
of Fig.lc and Id, the ratio of equilibrium constants of 


GR binding to the two chromatin states is a = Ka/K'^, 
and the ratio of chomatin opening and closing rates is 
P = ko/kc- 7 is not an equilibrium constant and is mod¬ 
ifiable but the cellular contents in remodeling enzymes 
and ATP. An energy-independent conformational equi¬ 
librium between PA and P'A would lead again to molec¬ 
ular trapping (of A and of the chromatin-remodeling en¬ 
zyme) . The equivalence between ko and a Poissonian rate 
is a gross approximation since the transition PA —> P'A 
encloses many elementary events, including the recruit¬ 
ment of enzymes, of ATP, catalytic steps etc, which will 
not be detailed here. Let be x and y the binding po¬ 
tentials of A and B respectively which are, for the non¬ 
equilibrium mechanisms x = A'(j[A], y = Kb[B\ and for 
the equilibrium mechanisms x = Ka [A] and y = Kb [B] . 
Fractional activity can be defined with these symbols and 
used for drawing 3D plots. They are listed below for the 
different models. 

6.1 Independent binding of A and B to 
the promoter 

Eq.(l) yields 


F = 


1 

1 -I- 7 



1 + 2 / 


( 12 ) 


The corresponding curve is shown in Fig.2a, when 
saturating A can trigger only one quarter of maximal 
activation ( 7 = 3). 


6.2 Putative equilibrium hierarchical 
model of Fig.la (Eq.(3)) 

x(l + i'(l + 7)) 

(1 -I- 7)(1 -I- a; -I- xy) 

The corresponding plot is shown in Fig.2b for 7 = 3. 


6.3 Putative equilibrium hierarchical 
model of Fig.lb (Eq.(4)). 

x{l P y) +-ty{a P x) 

(1 -I- 7)(1 -I- a; -I- y{a + x)) 

6.4 Non-equilibrium hierarchical model 
of Fig.lc (Eq.(ll)). 

In this system, chromatin relaxation can occur during 
the periods of absence of B, irrespective of whether A is 
present or not. 


F = 


ax (i+/3a;(l^) 
(1 + 7)[1 -I- aa;(l ■ 


Piy 
PA + y))] 


(15) 
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Figure 2. Comparative shapes of promoter activity curves in linear coordinates (left panels) and in Hill coordinates (right 
panels), x and y are the binding potentials of the TFs A and B used in the main text, and X and Y are their logarithms. The 
small inserts show 2D sections of the Hill plots at the indicated planes. In all cases, B is considered 3-time more potent than 
A for activating transcription, (a) The two TFs bind independently to the gene promoter (Eq.(12)). (b) Putative equilibrium 
hierarchical cooperativity model (Eq.(13)). (c) Hierarchical model with chromatin remodeling, in which the basal chromatin 
organization state can be restored only from the TF-free promoter (Eq.(18)), with the combination of parameters (a, P, 7 ) 
= (0.001, 2, 3). A and B are assumed to participate to the recruitment of transcription machineries in an additive manner. 


6.5 Non-equilibrium hierarchical model same development as in section 5, gives: 
of Fig.Id. 


ko[P] piA\P) = k4P']{l - piA\P')){l - piB\P')) (16) 


In this alternative possibility, chromatin closing to B can 
occur only for a TF-free promoter (specifically not from 
the P'A state). This possibility could for example be ex¬ 
plained by the persistent molecular association between 
A and chromatin-remodeling enzymes. In this case, the 


leading to the following steady state P/P' ratio 


^_ k,{l + KA[A]) _ 

[P'] koKA[A]{l + K'^[A]){l + KB[B]) 
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and to the fractional activity 

ax[l + I3x{l + y) + l3-fy{l + x)\ 

(l + ^)[l + a^(l + ;3(l + ^)(l + y))] ^ ^ 

A representative plot of this condition is shown in 
Fig.2c for (a, /3, 7 ) = (0.05, 2, 3). This set of param¬ 
eters is chosen to agree with experimental observations. 
Indeed, B is considered as more potent than A because 
it is in fact not a single TF, but a combination of sev¬ 
eral potent TFs (NFl and Oct). The higher affinity of 
A (GR) for the promoter (K^ > Ka), is suggested by 
the fact that the chromatin configuration permissive to 
NFl/Oct binding, strongly favours GR binding (Belikov 
et ah, 2004). 

7 Functional opportunities of¬ 
fered by the system 

7.1 A source of nonlinearity for structur¬ 
ing gene networks 

Sigmoidal genetic responses are the typical ingredients of 
multistable dynamic gene regularory networks. In a sim¬ 
ple example, if a gene subject to this mode of transcrip¬ 
tional regulation stimulates its own expression through a 
positive feedback, then, the sigmoidal and saturable ex¬ 
pression curve crosses twice the non-saturable first-order 
degradation line, thus generating bistability (Gherry and 
Adler, 2000). This general role of nonlinear responses in 
the formation of Boolean-like networks is not specific to 
hierarchical cooperativity and will not be detailed further 
here. 

7.2 Cooperativity between non¬ 
interacting TFs 

The mechanism of transcriptional cooperativity most 
widely reported and modeled in the literature is that 
obtained with TFs capable of physically interacting with 
each other and binding to a series of non-consensual DNA 
elements in a promoter. By this way, the direct inter¬ 
actions between the TFs help them to bind together to 
DNA, whereas their individual affinity for their DNA ele¬ 
ments would not have allowed their independent binding. 
But a recent study suggested that this mode of cooper¬ 
ativity is in fact doubtful (Ghu et ah, 2009), since di¬ 
rect interactions between adjacent TFs would lead to the 
clustering of TFs on DNA which inevitably contains non¬ 
specific TF binding sites. Indeed, all the TFs have a min¬ 
imal non-specific affinity for DNA (at least electrostatic). 
This point is interesting since it suggests that many in¬ 
teractions experimentally shown between TFs to explain 
cooperativity, could result from experimental drawbacks 
in detecting protein interactions (Mackay et ah, 2007). 


This problem no longer holds for the model of hierarchi¬ 
cal cooperativity described here. Moreover, the number 
of chromatin-remodeling machineries in the cell (Rippe 
et ah, 2007), which is so far intriguing, further supports 
the general importance of the present proposal. 

7.3 Participation of constitutively ex¬ 
pressed TFs to conditional expres¬ 
sion 

The NFl and Oct-1 TFs are generally expressed at high 
level by the laboratory cell lines used in the MMTV ex¬ 
periments cited above. However, in spite of their con¬ 
stitutive presence, the expression of MMTV integrated 
in ordered chromatin, is very low in absence of gluco¬ 
corticoid hormone (personal data not shown). Hence, 
NFl and Oct-1 contribute to the MMTV transcriptional 
strength but not to the decision to transcribe or not to 
transcribe. 

7.4 Adjustment of the degree of sensi¬ 
tivity and responsiveness of the pro¬ 
moter 

Expectedly, when chromatin remodeling is inhibited 
(/3 = 0), Eqs.(15) and (18) reduce to a simple hyper¬ 
bola Q;a ;/(1 -|- ax), but with chromatin remodeling, the 
behaviour of this system is unusual compared to classical 
modes of cooperativity. It can generate an ’’interrupter¬ 
like” mode of promoter functioning with strong non¬ 
linearity. For non-zero a and (3x, Eq.(18) can approach 
the Hill-like equation a/3(l -|- j)x‘^y/{l + q;/ 3(1 -f j)x‘^y), 
where the square exponent of x in absence of any pos¬ 
tulated dimerisation, reflects its dual participation in 
regulating both the P and P' promoter states. Interest¬ 
ingly, the self-cooperativity of x can be obtained even 
when B is transcriptionally inactive (7 = 0 ). In classical 
equilibrium mechanisms of promoter occupation, the re¬ 
sponsiveness and sensitivity to increasing concentrations 
of activated TFs, is fixed by the physico-chemical cellular 
conditions and by the values of equilibrium constants, 
which are themselves dictated by macromolecule struc¬ 
tures. For example in the case of TF dimerisation, the 
degree of cooperativity is not modifiable and determined 
by the affinity constants for given DNA elements. By 
contrast, in the present system, further adjustments are 
possible by tuning the activity and quantity of SWI/SNF 
chromatin remodeling enzymes (Fig.3) and the energy 
status of the cell (ATP), which influence altogether the 
j3 parameter. In particular, at given TF binding poten¬ 
tials, the chromatin-remodeling rate provides a precise 
control lever of the degree of promoter activation which 
can be coupled to additional threshold effects such as 
the buffering of low levels of mRNA expression by small 
RNAs, as shown in bacteria (Levine and Hwa, 2008). 
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Figure 3. Alternative view of the 3D plot of Fig.2c, showing 
the sigmoidal response of MMTV expression to the couple of 
activated TFs GR and Oct. For a given level of activity of 
SWI/SNF, this shape of the transcriptional activity surface 
makes the response of MMTV to GR and Oct strongly over¬ 
additive. Its near-horizontal slope at low GR and Oct values, 
allows to avoid inappropriate activation. Moreover, MMTV 
expression can be further adjusted by regulating the amount 
or the activity of the SWI/SNF enzymes and by other thresh¬ 
old effects which can cancel transcription activation by GR 
alone. 

Thought they are functionally important, the sub¬ 
tle differences between sigmoidal responses cannot be 
easily evaluated by eye in linear coordinates. Among 
the different mathematical tools developed for analysing 
non-hyperbolic fractional curves, the Hill representation 
long proved very useful because it allows to focus on 
the specificities of the systems (Cornish-Bowden and 
Koshland, 1975; Dahlquist, 1978). Using the logarithm 
of TF binding potentials allows to give to the range of 
binding potential between 0 to 1 , the same importance 
than between 1 and infinite, for better visualizing the ef¬ 
fects of ligand concentrations below midsaturation (for 
ln(Ar[TF]) = 0). The logarithm of the ratio of frac¬ 
tional activity vs inactivity ln(F/(l — F)) (’’logit” coor¬ 
dinate), allows to finely appreciate the behaviour of the 
system, including the degree of cooperativity between the 
TF(s) through the slope of the curve. An hyperbolic phe¬ 
nomenon gives a slope of 1. Indeed, when F = xj (1 -bee), 
then F/{1 — F) = x, thus eliminating saturation effects. 
In addition, in multidimensional Hill plots (3D in the 
present case), the relative participation of the different 
actors in the course of saturation can be visualized. Al¬ 
though the Hill plots are generally used for equilibrium 
phenomena such as hemoglobin oxygenation, they can 
also apply to steady states. The Hill equations corre¬ 
sponding to Eqs.(15) and (18) are Eqs.(19) and (20) re¬ 
spectively: 


ae 


X 


H, 


{X,Y) 


= In ■ 


1 + /3e^ (i^ei) + 


1 -f 7 -b ae 
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l + p[\^)+h 


(19) 


_ ae^[l-b/3e^(l-b e"*^)-b/37e^(l-b e^)] 
(X-Y) 2 -b 7 -b Q!e ^[7 -b /3(1 + e^) -b + e^)] 

( 20 ) 


where X = ln(a:) and Y = ln(?/). The correspond¬ 
ing plots are shown in the right panels of Fig.2 using 
the parameter combination (a, (3, 7 ) = (10“^, 2, 3). 
In these Hill surfaces, slopes of 1 correspond to free 
random (hyperbolic) binding, while non-unity slopes de¬ 
note the existence of collective influences in the system. 
Specifically, steep slopes reflect a phenomenon of coop¬ 
erativity increasing the sensitivity of the system to slight 
changes in ligand concentration. Near horizontal slopes 
and plateaus indicate the regions of relative TF ineffi¬ 
cacy as long as the concentration of the other TF is lim¬ 
iting. This is a situation of negative cooperativity. These 
Hill landscapes highlight the differences between the ba¬ 
sic (but doubtful) hierarchical mechanism (right panel 
of Fig.2b) and the nonequilibrium model (right panel of 
Fig.2c). While there is no limitation other than satura¬ 
tion in the response to large Y in the equilibrium model, 
this is no longer the case when Y > X in Fig.2c. In this 
respect, the latter model recovers some features of the 
independent system in which parallel increases of X and 
Y are necessary to allow their action. This property is re¬ 
lated to the fact that A can always escape the promoter 
and is not trapped contrary to the equilibrium model. 
This difference could be used as a tool for experimentally 
probing hierarchical systems. The hierarchical nature of 
these systems is illustrated by the preponderant role of 
A at low fractional activity. Hence, the active chromatin 
remodeling mechanism described here allows pronounced 
non-linearity, even for monomeric TFs, which can be fur¬ 
ther enhanced by other modes of cooperativity. 


7.5 Cumulating the layers of cooperativ¬ 
ity 

Hierarchical cooperativity provides an exquisite mode 
of sigmoidicity, in equilibrium (Fig.la) as well as non¬ 
equilibrium conditions (Fig.lc,d). In the equilibrium sys¬ 
tem, joint sigmoidicity is obtained only in the A + B 
bisector, by intersecting two series of orthogonal hyper¬ 
bolas (Fig.2b, illustrated by the 2D Hill curve at A = Y). 
In addition, in the active remodeling model, the response 
to A alone is also sigmoidal (visible along the A axis in 
the 3D plot of Fig.2c). The self-cooperativity of A fur¬ 
ther enhances the steepness of the global response in the 
bisector [A + B) (small 2D plot in Fig.2c). The max¬ 
imal Hill coefficients {tih) for the different models are, 
for the independent TFs of Fig.2a (Eq.(ll)): nniA) = 1, 
nuiB) = 1, uh^A + B) = 1; for the equilibrium model 
of Fig.la: n//(A) = 1, n_y(H) = 1, nniA B) = 2 
and for the non-equilibrium model of Fig.Id and Fig.2c: 
nniA) = 2, nniB) = 1, nniA + B) = 3. The sigmoidic¬ 
ity of this latter situation is illustrated in Fig.3. This 
source of sigmoidicity can surimpose to other ones, in¬ 
cluding: i) TF multimerization (neglected here) and ii) 
the cooperative recruitment of transcription machineries 
by DNA-bound TFs (Michel, 2010), rarely considered in 
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transcription modeling studies. For simplicity, it has not 
been not taken into account in the present study and 
Eq.(l) describes additive contributions of the TFs A and 
B to the global promoter activity. 


8 Speculative role of the present 
mechanism for the MMTV 

Certain mouse strains contain without apparent trouble, 
genome-integrated MMTV which are vertically trans¬ 
mitted over generations. MMTV-infected mouse cells 
can also remain healthy. Though viruses are generally 
detrimental for infected cells, the issue of an infection 
for the host cells often depends on the conditions. As 
it is counter-productive for a stowaway to destroy his 
vehicle, viral infection is not necessarily lytic. Indeed, 
host genome-integrated viruses have the opportunity to 
propagate passively as furtive aliens, through the mere 
spreading of the host cells. Accordingly, they developed 
strategies during evolution to preserve host viability as 
long as living conditions are satisfactory. In turn, when 
the viability of the host cells is menaced, the lytic phase 
is triggered and leads to the production of metabolically 
inert viral particles which are more resistant to delete¬ 
rious conditions. This strategy has been observed in 
the prokaryotic world, for example in the case of the 
lambda bacteriophage in lysogenic bacteria, but it can 
also apply to certain eukaryotic integration viruses, such 
as the MMTV provirus which generally remains dor¬ 
mant in adults, unless they are submitted to stresses. 
MMTV expression is weak in stressless conditions since 
nucleosomes ensure its transcriptional silencing. In cells 
containing GR, glucocorticoid hormones trigger MMTV 
expression. Glucocorticoid hormones (corticol, corticos¬ 
terone), are the hormones of nervous stress, which acti¬ 
vate the whole panoply of GR activities (nuclear import, 
DNA binding, transactivation, recruitement of BRGl). 
The secondary TFs which strongly enforce the GR ac¬ 
tion are NFl but also Oct-1 or Oct-2. Interestingly, 
the preferential binding site for Oct-2 defined in (Rhee 
et ah, 2001) precisely corresponds to the Oct module 
present in the MMTV promoter. It is inducible by bacte¬ 
rial lipopolysaccharides (LPS) and inflammatory signals. 
Hence, several types of stresses: nervous (GR) and infec¬ 
tious (Oct-2), concur to activate MMTV. The sigmoidal 
shape of the response shown in Fig.2c, is such that the 
combination of the two types of stress is required to trig¬ 
ger transcription. Moreover, the near horizontal slope 
of the transcription surface and the very low responsive¬ 
ness to low GR and Oct concentrations (Fig.3), allow to 
buffer stochastic fluctuations of these TFs. By this way, 
MMTV can remain latent in moderately stressed cells, 
and is revived upon conjunction of stresses (Fig.4). 


GRE Oct 



Figure 4. Example of hierarchical transcriptional cooper- 
ativity mediated by chromatin remodelling, (a) Schematic 
representation of the proximal MMTV promoter critically 
regulated by a nucleosome (grey), which is positioned to 
overlap the binding sites for GR and Oct-2. In this configu¬ 
ration, only GR can bind its target site, owing to its capacity 
to interact with nucleosomal DNA. Its fixation then triggers 
the recruitment of the chromatin remodeller BRGl, which in 
turns allows the fixation of the Oct-2 factor requiring a fully 
accessible DNA helix, (b) This systems predicts an overad¬ 
ditive combination effect and the response of the MMTV to 
multiple stresses. 

This could be the case for example when the host¬ 
ing mouse is both frightened, with production of gluco¬ 
corticoid hormone (for example if a cat appears in the 
neighbourhood) and wounded (leading to a bacterial in¬ 
fection and to Oct-2 induction). When these conditions 
are reunited, the mouse’s life is probabilistically compro¬ 
mised and it is beneficial for the MMTV to escape it 
before sinking with it. MMTV expression can be partic¬ 
ularly important in lymphocytes because these cells are 
cellular reservoirs for MMTV, undergo apoptosis upon 
glucocorticoid exposure and display strong Oct-2 induc¬ 
tion by inflammatory stress (Bendall et ah, 1997). This 
transcriptional arbitration is equivalent to that of a crisis 
board but is more economic. The conversion of random 
interactions into discerning actions is a typical character¬ 
istic of dissipating systems, involving in the present case 
energy-dependent chromatin remodeling. 
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9 Conclusion 

The model proposed here is a simplification omitting 
many actors in the MMTV promoter story, but is suffi¬ 
cient to reconcile conflicting data. While the first articles 
convincingly demonstrated that GR binding opens the 
way to NFl and Oct-1, further studies showed that NFl 
and Oct-1 are present prior to glucocorticoid hormone 
addition (Belikov et ah, 2004). The presetting action 
of NFl suggested in this latter article was interpreted 
as a locking action of NFl, that was suggested to click 
nucleosome positioning in a unique configuration. This 
interpretation is fully consistent with the mechanism pro¬ 
posed here, in which chromatin closing and NFl binding 
are mutually exclusive events. In the present model, a 
fuzzy pattern is expected if the P and P' states and their 
transition intermediates coexist in the cell population. 
This coexistence is possible in presence of GR alone (lig¬ 
and A in Fig.lc), but not of NFl or Oct-1 (ligand B in 
Fig.lc). The P state corresponds to the positioning of 
nucleosomes thermodynamically favoured by nucleotide 
sequence-specific DNA bendability (Pina et ah, 1990). 
The P' state is a less stable configuration, whose forma¬ 
tion is forced by SWI/SNF ATPases and which is then 
locked by NFl/Oct-1 as long as present. This scheme is 
satisfactory in that it allows to explain previous obser¬ 
vations seemingly contradictory: (i) the initiation role of 
activated GR on NFl/Oct-1 fixation, (ii) the presetting 
action of NFl/Oct-1 on GR exchanges, (iii) the fact that 
GR is not trapped in presence of NFl/Oct-1. Consid¬ 
ering the abundance of chromatin remodeling factors in 
the cell (Rippe et ah, 2007), such a mechanism could 
be very general and provide a widespread mode of co- 
operativity between TFs that do not directly interact 
with each other. In addition, these enzymes render cell- 
context specific, the role of ubiquitous actors such as 
the DNA-binding elements for TFs that are common to 
several cell types. Two nuclear receptors: GR and PR 
(progesterone receptor), are of equivalent strength and 
share the same DNA modules in the MMTV promoter; 
but interestingly, in a cellular context permissive for 
GR, PR fails to activate MMTV integrated into ordered 
chromatin, but induces MMTV when transfected in an 
open chromatin state (Smith et ah, 1997). Accordingly, 
PR is unable to induce chromatin remodeling at stably 
integrated MMTV templates in these cells (Smith et 
ah, 1997; Fryer and Archer, 1998) and the reciprocal 
situation is obtained in other cellular contexts (T47D, 
personal data). The mechanism proposed in this study 
could be a pivotal device for the management of the 
eukaryotic genomes based on their nucleosomal organ¬ 
isation. It allows: (i) to solve apparent discrepancies 
between experimental observations, so far barely rec¬ 
oncilable in equilibrium conditions; (ii) to establish a 
primary and highly tunable mode of cooperativity be¬ 
tween TFs, considering that the chromatin-remodeling 
enzymes are themselves subject to refined regulations; 


(iii) and to bypass the need for direct interactions be¬ 
tween them, which is questioned in (Chu et ah, 2009). 

Revision of the concept of pioneer transcrip¬ 
tion factors. Pioneer transcription factors are defined 
as developmental factors opening the way to secondary 
transcription factors. In this sequential mode of ac¬ 
tion, the pioneer transcription factor is envisioned as au¬ 
tonomous whereas the secondary transcription factor is 
tributary of the pioneer one. For example, the pioneer 
factors FOXAl (also involved in the MMTV system), 
AP 27 , PBXl and GATA3, are supposed to preset chro¬ 
matin and allow the fixation of the estrogen receptor-a 
(ERa) in mammary lumenal epithelial cells. More than 
80% of the ERa-binding sites are associated to the fixa¬ 
tion of one of these pioneer factors (Magnani and Lupien, 
2014). The depletion of these factors prevents ERa from 
binding, but the reverse has recently also been shown 
true (Caizzi et ah, 2014), which singularly challenges the 
unilateral dependence of secondary transcription factors 
on pioneer factors. By contrast, the present model is 
compatible with a reciprocal dependence between these 
factors. As shown in Fig.la,b, the secondary factor can 
strengthen the fixation of the pioneer factor by extending 
the fraction of time of the remodeled chromatin state, to 
which the pioneer factor can bind with a higher affinity. 
In addition, this mechanism is dynamic, contrary to a 
static hierarchical view FOXAl —)■ ERa, which is poorly 
compatible with the half residence time of FOXAl of 4 
minutes, as measured by fluorescence microscopy (Sekiya 
et ah, 2009). The mechanism described here predicts mu¬ 
tual influences between transcription factors, creating a 
key combination effect for turning on or off gene expres¬ 
sion. 

The present manuscript is an extended version of the 
article: Hierarchical cooperativity mediated hy chromatin 
remodeling; the model of the MMTV transcription regu¬ 
lation. Michel, D. 2011. J. Theor. Biol. 287, 74-81. 
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